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Abstract 

Roe and Woodroofe (RW) have suggested that certain conditional probabil- 
ities be incorporated into the "unified approach" for constructing confidence 
intervals, previously described by Feldman and Cousins (FC). RW illustrated 
this conditioning technique using one of the two prototype problems in the 
FC paper, that of Poisson processes with background. The main effect was 
on the upper curve in the confidence belt. In this paper, we attempt to apply 
this style of conditioning to the other prototype problem, that of Gaussian 
errors with a bounded physical region. We find that the lower curve on the 
confidence belt is also moved significantly, in an undesirable manner. 

PACS numbers: 06.20.Dk, 14.60.Pq 
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I. INTRODUCTION 



Roe and Woodroofe |jj have made an interesting suggestion for modifying the "uni- 
fied approach" to classical confidence intervals which Feldman and I advocated in Ref. [Q. 
They invoke the use of "conditioning", namely replacing frequentist coverage probabilities 
with conditional probabilities, still calculated in a frequentist manner, but conditioned on 
knowledge gained from the result of the particular experiment at hand. 

Roe and Woodroofe (RW) illustrate their suggestion using one of the two prototype prob- 
lems, that of Poisson processes with background. Suppose, for example that an experiment 
observes 3 events (signal plus background). Then the experimenters know that, in that 
particular experiment, there were 3 or fewer background events. RW therefore calculate the 
frequentist coverage using an ensemble of experiments with 3 or fewer background events, 
rather than the larger unrestricted ensemble which we used. Thus, the RW ensemble changes 
from experiment to experiment. Conditioning on an equality has a long history in classical 
statistics. (Ref. [|I|] contains key references.) However, conditioning on an inequality, as RW 
do when the number of events is greater than zero, is perhaps less well founded, and it is 
interesting to explore the consequences. 

In this paper, we attempt to apply RW-like conditioning to the other prototype problem, 
that of Gaussian errors with a bounded physical region. The result is similar to the Poisson 
problem analyzed by RW, but difficulties which were apparently masked by the discrete 
nature of the Poisson problem now arise. In particular, the lower endpoints of confidence 
intervals are moved significantly in an undesirable direction. 



II. THE UNIFIED APPROACH TO THE GAUSSIAN- WITH-BOUNDARY 

PROBLEM 

As in Ref. ||, we consider an observable x which is the measured value of parameter /i 
in an experiment with a Gaussian resolution function with known fixed rms deviation a, set 
here to unity. I.e., 

P(x\n) = -Lexp(-(;r - fi) 2 /2). (2.1) 

We consider the interesting case where only non-negative values for n are physically allowed 
(for example, if fi is a mass). 

The confidence-belt construction in Ref. @ proceeded as follows. For a particular x, 
we let /^best be the physically allowed value of [i for which P(x\fi) is maximum. Then 

/Ubcst = max(0,x), and 

n^Kcst) = { Vv f \ /9W nr (2-2) 
[ exp(— x /2)/v2vr, x < 0. 

We then compute the likelihood ratio R, 

= Pjjjj = f exp(-(x - ^) 2 /2), x > 
P(x\n heBt ) \ exp(x/i - fi 2 /2), x<0. 
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During our Neyman construction of confidence intervals, R determines the order in which 
values of x are added to the acceptance region at a particular value of /x. In practice, this 
means that for a given value of /i, one finds the interval [xi,^] such that R(xi) = R{x2) 
and 

/ P(x\fi)dx = a, (2.4) 

J xi 

where a is the confidence level (C.L.). We solve for x\ and X2 numerically to the desired 
precision, for each \i in a fine grid. With the acceptance regions all constructed, we then 
read off the confidence intervals [/ii,/^] as in Ref. 0. 



III. INVOKING CONDITIONING IN THE GAUSSIAN- WITH-BOUNDARY 

PROBLEM 

In order to formulate the conditioning, we find it helpful to think of the measured value 
x as being the sum of two parts, the true mean fi t and the random "noise" which we call e: 

x = /j, t + e. (3.1) 

We are considering the case where it is known on physical grounds that /it > 0. Thus, if an 
experimenter obtains the value Xq in an particular experiment, then he or she knows that, 
in that particular experiment, 

s < x . (3.2) 

For example, if the experimenter measures /j and obtains xq = —2, then the experimenter 
knows that e < —2 in that particular experiment. This information is analogous to the 
information in the Poisson problem above in which one knows that in the particular experi- 
ment, the number of background events is 3 or fewer. We thus use it the manner analogous 
to that of RW: our particular experimenter will consider the ensemble of experiments with 
e < xq when constructing the confidence belt relevant to his or her experiment. 

We let P(x\fj,,e < Xq) be the (normalized) conditional probability for obtaining x, given 
that e < Xq. In notation similar to that of RW, this can be denoted as q^°(x): 

g-(*)-P(^<a; ) = { (3.3) 

10, X ^> fl -\- Xq. 

Given xq, at each x we find //best, that value of \x which maximizes P(x\fi,e < xq): 

{x, xq > and x > 

x — Xq, Xq < and X > Xq (3.4) 
0, otherwise 

In the notation of Ref. M, P(x\/j l \ )e8t , e < x ) is then 



1, Xq > and x > 



max q x ?(x) = — x <^ exp(-x§)/2, x < and x > x (3.5) 

v V 27r(erf (V V2) + 1) [ exp (_ x 2 ) /2; otherwise 
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Then the ratio R of Eqn. [2.3| is replaced by 



R xo (li,x) 




(3.6) 



iriax^/ q x ? (x) ' 



which vanishes if x > fi + Xq, and otherwise is given by 



R x °{^x) 



exp(— (x — /i) 2 /2), xq > and x > 

exp((— (x — yu) 2 + Xq)/2), xq < and x > x 
exp(x/i — /i 2 /2), otherwise 



(3.7) 



Figures [i] through |3] show graphs of q x °(x), max^/ q x ?(x), and R x °(fi,x), for three values of 
y^, for each of three values of x . 

We let Cj; (/x) be the value of c for which 



x:R x 0(p,x)<c 




The modified confidence interval consists of those /i for which 



R x °(fi,xo) > c xo (n). 



(3.9) 



Note that this entire construction depends on the value of xo obtained by the particular 
experiment. An experiment obtaining a different value of xq will have a different function 
in Eqn. |3^, and hence a different confidence belt construction. Figure |] shows examples of 
such constructions for six values of Xq. The vertical axis gives the endpoints of the confi- 
dence intervals. Each different confidence belt construction is used only for an experiment 
obtaining the value xq which was used to construct the belt. The interval [/ii,/i2] at x = xo 
is read off for that experiment; the rest of that plot is not used. 

Finally, we can form the graph shown in Fig. [| by taking the modified confidence interval 
for each x$, and plotting them all on one plot. These are tabulated in Table |, which includes 
for comparison the unconditioned intervals from Table X of Ref . . 

Fig. |6] shows the modified intervals plotted together with the unified intervals of Ref. [0 . 
The modified upper curve is shifted upward for negative x, which results in a less stringent 
upper limit when e is known to be negative; this feature is considered desirable by some. 
The lower curve, however, is also shifted upward: for all xq > 0, the interval is two-sided. 
We find this to be a highly undesirable side-effect. 

It is interesting to consider what happens if one applies Fig. ^ to an unconditioned 
ensemble. The result can be seen by drawing a horizontal line at any /i in Fig. || and inte- 
grating P(x\fi) (Eqn. pTTD along that line between the belts. For small /i, there is significant 
under cover age, while for /z near 1.0, there is significant overcoverage. The undercoverage 
was surprising, since the conditioned intervals always cover within the relevant subset of the 
ensemble. However, conditioning on an inequality means that these subsets are not disjoint. 

The undesirable raising of the lower curve is present in the Poisson case, as can be seen 
in Figure 1 of Ref. ||]. However, there the discreteness of the Poisson problem apparently 
prevents the curve from being shifted so dramatically, and the two-sided intervals do not 
extend to such low values of the measured n. 
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IV. CONCLUSION 

In this paper, we apply conditioning in the style Roe and Woodroofe to the Gaussian- 
with-boundary problem. We find that the transition from one-sided intervals to two-sided 
intervals undesirably moves to the origin. This reflects a general feature of confidence interval 
construction: when moving one of the two curves, the other curve moves also. In the Poisson- 
with-background problem, the undesirable movement was not large, but in the Gaussian- 
with-boundary problem, the effect is quite substantial. 
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FIG. 1. Graphs of q^°(x) (top row), max^ q*?(x) (middle row), and R x °(fi,x) (bottom row), 
for ix = 0.01. The columns are for xo = —1, 0, and 1. Each graph in the bottom row is the quotient 
of the two graphs above it. 
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FIG. 2. Graphs of q^°(x) (top row), max^ qp{x) (middle row), and R x °(/j,,x) (bottom row), 
for \x = 0.5. The columns are for xq = —1, 0, and 1. Each graph in the bottom row is the quotient 
of the two graphs above it. 
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FIG. 3. Graphs of q^°(x) (top row), max^/ qp(x) (middle row), and R x °(/j,,x) (bottom row), 
for fj, = 2.5. The columns are for xq = —1, 0. and 1. Each graph in the bottom row is the quotient 
of the two graphs above it. 
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FIG. 4. Conditional confidence belts for the six sample values of x$ indicated. Each plot is 
used only to find the [/ii,^] interval at x equal to the xq used to construct it; that interval is 
transferred to Figj5[ 
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FIG. 5. Plot of RW-inspired 90% conditional confidence intervals for mean of a Gaussian, 
constrained to be non-negative, described in the text. 
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FIG. 6. Plot of RW-inspired 90% conditional confidence intervals (solid curves) , superimposed 
on the unconditioned intervals of Ref. 0] (dotted curves). 
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TABLES 



TABLE I. 90% C.L. confidence intervals for the mean fi of a Gaussian, constrained to be 
non- negative, as a function of the measured mean xq, for the RW conditioning method, and for 
the unified approach of Feldman and Cousins. All numbers are in units of a. The conditioned 
numbers may be inaccurate at the level of ±0.01 due to the computational grid used. 
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